The Color Distributions of Globular Clusters in Virgo Elliptical Galaxies

This Letter presents the color distributions of the globular cluster (GC) systems of 12 Virgo elliptical galaxies, measured using data from the Hubble Space Telescope. Bright galaxies with large numbers of detected GC's show two distinct cluster populations with mean V-I colors near 1.01 and 1.26. The GC population of M86 is a clear exception; its color distribution shows a single sharp peak near V-I=1.03. The absence of the red population in this galaxy, and the consistency of the peak colors in the others, may be indications of the origins of the two populations found in most bright elliptical galaxies.Comment: 5 pages, 1 figure, to be published in ApJ Letters Corrections to introductio

Massively parallel simulations of binary black holes with Dendro-GR

We present results from the new Dendro-GR code. These include simulations of binary black hole mergers for mass ratios up to q=16. Dendro-GR uses Wavelet Adaptive Multi-Resolution (WAMR) to generate an unstructured grid adapted to the spacetime geometry together with an octree based data structure. We demonstrate good scaling, improved convergence properties and efficient use of computational resources. We validate the code with comparisons to LazEv

The Surface Brightness Fluctuations and Globular Cluster Populations of M87 and its Companions

Using the surface brightness fluctuations in HST WFPC-2 images, we determine that M87, NGC 4486B, and NGC 4478 are all at a distance of ~16 Mpc, while NGC 4476 lies in the background at ~21 Mpc. We also examine the globular clusters of M87 using archived HST fields. We detect the bimodal color distribution, and find that the amplitude of the red peak relative to the blue peak is greatest near the center. This feature is in good agreement with the merger model of elliptical galaxy formation, where some of the clusters originated in progenitor galaxies while other formed during mergers.Comment: 5 pages, 2 figure

Relativistic MHD with Adaptive Mesh Refinement

This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference Convex ENO method (CENO) in 3+1 dimensions, and the AMR is Berger-Oliger. Hyperbolic divergence cleaning is used to control the $\nabla\cdot {\bf B}=0$ constraint. We present results from three flat space tests, and examine the accretion of a fluid onto a Schwarzschild black hole, reproducing the Michel solution. The AMR simulations substantially improve performance while reproducing the resolution equivalent unigrid simulation results. Finally, we discuss strong scaling results for parallel unigrid and AMR runs.Comment: 24 pages, 14 figures, 3 table

Perturbed disks get shocked. Binary black hole merger effects on accretion disks

The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk's angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk's angular momentum, the disk's dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.Comment: 10 pages, 13 figures. v2: Minor changes, version to be published in PR

High-Fidelity Multidisciplinary Design Optimization Methodology with Application to Rotor Blades

A multidisciplinary design optimization procedure has been developed and applied to rotorcraft simulations involving tightly-coupled, high-fidelity computational fluid dynamics and comprehensive analysis. A discretely-consistent, adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows on unstructured, dynamic, overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system for high-fidelity rotorcraft analysis is verified; simulation results exhibit good agreement with established solutions. A constrained gradient-based design optimization for a HART-II rotorcraft configuration is demonstrated. The computational cost for individual components of the multidisciplinary sensitivity analysis is assessed and improved

Simulating binary neutron stars: dynamics and gravitational waves

We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the \had computational infrastructure for distributed adaptive mesh refinement.Comment: 14 pages, 16 figures. Added one figure from previous version; corrected typo

Discovery of New Ultracool White Dwarfs in the Sloan Digital Sky Survey

We report the discovery of five very cool white dwarfs in the Sloan Digital Sky Survey (SDSS). Four are ultracool, exhibiting strong collision induced absorption (CIA) from molecular hydrogen and are similar in color to the three previously known coolest white dwarfs, SDSS J1337+00, LHS 3250 and LHS 1402. The fifth, an ultracool white dwarf candidate, shows milder CIA flux suppression and has a color and spectral shape similar to WD 0346+246. All five new white dwarfs are faint (g > 18.9) and have significant proper motions. One of the new ultracool white dwarfs, SDSS J0947, appears to be in a binary system with a slightly warmer (T_{eff} ~ 5000K) white dwarf companion.Comment: 15 pages, 3 figures, submitted to ApJL. Higher resolution versions of finding charts are available at http://astro.uchicago.edu/~gates/findingchart

The Dark Energy Survey Data Management System

The Dark Energy Survey collaboration will study cosmic acceleration with a 5000 deg2 griZY survey in the southern sky over 525 nights from 2011-2016. The DES data management (DESDM) system will be used to process and archive these data and the resulting science ready data products. The DESDM system consists of an integrated archive, a processing framework, an ensemble of astronomy codes and a data access framework. We are developing the DESDM system for operation in the high performance computing (HPC) environments at NCSA and Fermilab. Operating the DESDM system in an HPC environment offers both speed and flexibility. We will employ it for our regular nightly processing needs, and for more compute-intensive tasks such as large scale image coaddition campaigns, extraction of weak lensing shear from the full survey dataset, and massive seasonal reprocessing of the DES data. Data products will be available to the Collaboration and later to the public through a virtual-observatory compatible web portal. Our approach leverages investments in publicly available HPC systems, greatly reducing hardware and maintenance costs to the project, which must deploy and maintain only the storage, database platforms and orchestration and web portal nodes that are specific to DESDM. In Fall 2007, we tested the current DESDM system on both simulated and real survey data. We used Teragrid to process 10 simulated DES nights (3TB of raw data), ingesting and calibrating approximately 250 million objects into the DES Archive database. We also used DESDM to process and calibrate over 50 nights of survey data acquired with the Mosaic2 camera. Comparison to truth tables in the case of the simulated data and internal crosschecks in the case of the real data indicate that astrometric and photometric data quality is excellent.Comment: To be published in the proceedings of the SPIE conference on Astronomical Instrumentation (held in Marseille in June 2008). This preprint is made available with the permission of SPIE. Further information together with preprint containing full quality images is available at http://desweb.cosmology.uiuc.edu/wik